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1 Introduction. 



This paper is a continuation of our previous investigation on the truncated 
matrix trigonometric moment problem by the operator approach in pQ. The 
truncated matrix trigonometric moment problem consists of finding a non- 
decreasing CjVxAr-valued function M(t) = (jTOfc,z)^z=o> t £ [0, 2ir], M(0) = 0, 
which is left-continuous in (0,27r], and such that 



i-2n 

/ e int dM(t) = S n , n = 0,l,...,d, 
J o 



where {;S n }^ =0 is a prescribed sequence of (N x N) complex matrices (mo- 
ments). Here Af£N and d £ Z + are fixed numbers. Set 



T d = (Si 



-j)i,j=0 
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where 



S, := S*_ 



ki 



-d,-d+ 1 



and {S'nj^—Q are from ([T]). It is well known that the following condition: 

T d >0, 



(3) 



is necessary and sufficient for the solvability of the moment problem ([T]) 
(e.g. [2J). The moment problem ([T]) is said to be determinate if it has a 
unique solution and indeterminate in the opposite case. We shall omit here 
an exposition on the history and recent results for the moment problem (prj. 
All that can be found in pQ. 

The aim of our present investigation is to derive a Nevanlinna-type formula 
for the truncated matrix trigonometric moment problem (TMTMP) in a 
general case. Namely, we shall only assume that d > 1, condition ([3]) is 
satisfied and the moment problem is indeterminate, i.e. it has more than 
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one solution. The coefficients of the corresponding matrix linear fractional 
transformation are explicitly expressed by the prescribed moments. Notice 
that in some situations (e.g. during a multiple application) the Nevanlina- 
type formula has the advantage that the numbers of rows and the numbers 
of columns of its coefficients are less or equal to N. Easy conditions for the 
determinacy of the moment problem in terms of the prescribed moments are 
given. 

Notations. As usual, we denote by R, C, N, Z, Z + , the sets of real 
numbers, complex numbers, positive integers, integers and non-negative in- 
tegers, respectively; D = {z £ C : \z\ < 1}. The set of all complex matrices 
of size (rn x n) we denote by C mX n) rn,n £ N. If M 6 C mX n 

then M 1 

denotes the transpose of M, and M* denotes the complex conjugate of M. 
The identity matrix from C nxn we denote by I n , ti£N. 

If H is a Hilbert space then (•,•)# an d || ■ \\h mean the scalar product 
and the norm in H, respectively. Indices may be omitted in obvious cases. 
By we denote the finite-dimensional Hilbert space of complex column 
vectors of size N with the usual scalar product (x, y)cN = Ylf^Q 1 x jVj> f° r 
x,y e C N , x = (x ,xi, . . .,x N -i) T , y = (y ,yi, . . . ,y N -i) T , Xj,yj G C. 
For a linear operator A in H, we denote by D(A) its domain, by R{A) 
its range, by KerA its null subspace (kernel), and A* means the adjoint 
operator if it exists. If A is invertible then A^ 1 means its inverse. A means 
the closure of the operator, if the operator is closable. If A is bounded then 
1 1 .A 1 1 denotes its norm. For a set M C H we denote by M the closure of M 
in the norm of H. If M has a finite number of elements, then its number of 
elements we denote by card(M). For an arbitrary set of elements {x n } ne / 
in H, we denote by Lin{x n } ne / the set of all linear combinations of elements 
x n , and span{x n } n6 / := Lin{x n } ng /. Here I is an arbitrary set of indices. 
By Eh we denote the identity operator in H, i.e. Ehx = x, x G H. In 
obvious cases we may omit the index H. If Hi is a subspace of H, then 
Pjlj = Pjj is an operator of the orthogonal projection on H\ in H. 



2 The determinacy of the TMTMP. A Nevanlinna- 
type formula for the TMTMP. 

Let the moment problem (fT]), with d > 1, be given and condition ([3]), with 
Td from ([2]), be satisfied. Let 

T d = (7n,m)i d m=!^~\ S k = (S k; s,i)^tZ} , -d<k<d, 
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where 7 n , m , Sk ;s ,i G C Observe that 

7fcAf+s,r7V+i = Sk-r;s,l, < k,r < d, < s, Z < N - 1. (4) 

We repeat here some constructions from pQ. Consider a complex linear vec- 
tor space ft, which elements are row vectors u = (uq, «i, «2, •••) Wd4-iW-i)> 
with «„ £ C, < n < (d+ l)iV — 1. Addition and multiplication by a scalar 
are defined for vectors in a usual way. Set 

S n = (5n,0,5 n ,l A,2, ••• A,(d+l)AT-l), < fl < (d + 1)N - 1, 

where <5 njr is Kronecker's delta. In .f) we define a linear functional B by the 
following relation: 

(d+l)JV-l 
B(u,w) = a n b r -f n!r , 

n,r=0 

where 

(d+l)JV-l (d+l)iV-l 
?2 = fln^n) ^ = 6 r e r , a n ,b r £ C. 

n=0 r=0 

The space f) with 1? form a quasi-Hilbert space ([3]). By the usual procedure 
of introducing of the classes of equivalence (see, e.g. j3j), we put two elements 
u, w from fj to the same class of equivalence denoted by [u] or [w], if B(u — 
w,u — w) = 0. The space of classes of equivalence is a (finite-dimensional) 
Hilbert space. Everywhere in what follows it is denoted by H. Set 

x n :=[e n ], 0< n< (d + l)N- 1. 

Then 

(x n ,x m ) H = 7n,m, < n,m < (d + l)iV - 1, (5) 

and span^}^ 1 ^ -1 = L M^n}„t!o * = H - Set L ^ := ^{^n}^ 1 - 
Consider the following operator: 

dN-l dN—1 
Ax = y~] a k x k+N , x = OfcXfc, a fc G C. (6) 

fc=0 fc=0 

By [H Theorem 3] all solutions of the moment problem ([1]) have the following 
form 

M{t) = (m kJ (t))^ , te[0M, (7) 
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where rrikj are obtained from the following relation: 

J Y^-Ft dm kA t ) = (i E H-C(A®<l> c T 1 x k ,x j ) H , (GD. (8) 

Here is an analytic in D operator-valued function which values are linear 
contractions from H Q D{A) into H QR(A). Conversely, each analytic in D 
operator-valued function with above properties generates by relations 0- 
(jHI a solution of the moment problem (pQ). The correspondence between 
all analytic in B operator-valued functions with above properties and all 
solutions of the moment problem fll]) is bijective. 

Since we are going to obtain a Nevanlinna-type formula for the indeter- 
minate TMTMP, it is important to obtain some easy necessary and sufficient 
conditions for the determinacy of the TMTMP. 

Theorem 1 Let the moment problem UP, with d > 1, be given and condi- 
tion with Td from (0|) ; be satisfied. Let the operator A in the Hilbert 
space H be constructed as in (G>|). The following conditions are equivalent: 

(A) The moment problem (OP is determinate; 

(B) The defect numbers of A are equal to zero; 

(C) For each fixed number r, dN < r < dN + N — 1, the following linear 
system of simultenuous equations: 

dN-l 

a r,n7n,j = lr,ji (9) 

n=0 

with unknowns a r o,avi) ■ ■ ■ i a rdN-li has a solution. Here the num- 
bers 7 V are defined by 

If the above conditions are satisfied then the unique solution of the mo- 
ment problem (CP is given by the following relation: 

M(t) = (m k ,j(t))%J= , m kJ {t) = (Etx k ,Xj) H , (10) 

where E t is the left- continuous orthogonal resolution of unity of the unitary 
operator A, which is piecewise constant. 

Proof. (A)=^(B). First, we notice that the defect numbers of A are 
always equal, because H is finite-dimensional and A is isometric. If the 
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defect numbers are greater then zero, then we can choose unit vectors u\ G 
H D(A) and u 2 G H Q R(A). We set §^cu\ = eu 2 , Vc G C, and §qu = 0, 
u e H Q D(A): u Linjui}; £ G D. On the other hand, we set $^ = 
0. Functions <E>£ and <3?£ produce different solutions of the TMTMP by 
relation ([8]). 

(B)=^(A). If the defect numbers are zero, then the only admissible function 
$ in relation §8§ is <3?^ = 0. 

Observe that (B) ^ (D{A) = H) & 



XdN, XdN+1, ■■■ , XdN+N- 



XdN = l^n=0 a dN,nX n 
^dN-1 
XdN+1 — Z^n=0 a dN+l,n x n 



.! G Lin{x n }^ " X ) 



\ 



OtdN.n, OidN+U 



,OidN+N-l,n S 



Z^n=0 a dN+N-l,nX n 



(x dN ,Xj) H = (En=0 1 a dN,nX n ,Xj) H 
(xdN+l,Xj)H = (Sn=0 "dW+l.n-^nj^Off 



V 



(.XdN+N-l>Xj)H = (Dn=0 1 a dN+N-l,nXn,Xj)H 
where OCdN,ni a dN+X,n: ■ ■ ■ ) a dN+N— l,n 

£ C, < j < dN + N - 1) ^ (C). □ 
We shall continue our considerations started before the statement of 
Theorem [TJ In what follows we assume that the TMTMP is indeterminate 
and the both defect numbers of A are equal to 5 = S(A), 5 > 1. 
Let us apply the Gram-Schmidt orthogonalization procedure to the vectors 
xq,xi, . . . ,XdN+N-i- During this procedure we shall use the numbers 7 V 
defined by (jl]) and the property ([5]). 
Step j; < j < dN + N - 1. Calculate 



^2 (xj,yk)HVk 

k: 0<fc<j-l, n fc ^0 



H 



where the sum on the right can be empty. If rij ^ 0, then we set 
1 



Vj 



^2 (.Xj,y k ) H y k 



(12) 



k: 0<Jfe<j-l, n k ^0 



If rij = 0, we pass to the next step. 
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Remark 1 Notice that there exists a nonzero nj with < j < N— 1. In the 
opposite case we would have \\xjWjj = =0, < j < N — 1. Therefore 
So = and M(t) = 0, an<i i/iis contradicts to the indeterminacy of the 
TMTMP. 

Remark 2 By ilty) every yj can be expressed as a linear combination of 
xq, x%, . . . , Xj. Thus, numbers nj can be calculated using the prescribed mo- 
ments by relations |2p and J^J). 

Set Oi = {j : < j < dN + N - 1, nj / 0}. Then 21 := {y^je^ 
is an orthonormal basis in H. Moreover 2li := {yj}j^n i: j<N-i is an or- 
thonormal basis in L^r, and 21-2 := {yj)j&Q,y. j<dN-i is an orthonormal basis 
in Lmixn}^ 1 = D(A). Therefore 2l 3 := {yj}jenv dN<j<dN+N-i is an 
orthonormal basis in H © D{A). Consequently, 5 < N. 
Observe that card(2la) = 5 > 1. Set p := card(2li), r := card(2l2). Notice 
that 1 < p< N, t> p>l. 

The k-th element, counting from zero, of the set 21, arranged in the order of 
construction of its elements, we denote by Uk, k = 0, 1, ...,r + 5 — 1. Then 

21 = Mitt 1 ' = WCJ, 2i 2 = 2i 3 = K-}^ 1 . 

We need one more orthonormal basis in H. Observe that 21^ := { v k}k=o> 
where Vk '■= Auk, is an orthonormal basis in R(A). Notice that R(A) = 
Linjxn}^^" 1 . Therefore the linear span of vectors {vkYk=o^ {xn}n=o> 1S 
equal to H. Hence, the linear span of vectors {vk} T k Z^, {u^^Zq, is equal to 
H , as well. 

Let us apply the Gram-Schmidt orthogonalization procedure to the vectors 
vo, v 1, . . . , w T _i, no, Mi, ... , u p -\. Like in the above procedure, we shall use 
the numbers 7 V defined by ([4]) and the property ([5]). Observe that the first 
r elements are already orthonormal. 
Step j; < j < p — 1. Calculate 

T-l 

Uj -^2(Uj,Vi) H Vl- ^ (Uj,fk)Hfk , (13) 

1=0 k: 0<k<j-l, m k ^0 

where the last sum on the right can be empty. If nij 7^ 0, then we set 

fj ■= — I u i ~ J2( u j, v j)H v j ~ Yl (uj,fk)Hfk ] • (14) 

J y 1=0 k: 0<k<j~l, m fe ^0 J 

If mj = 0, we pass to the next step. 



H 
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Set n 2 = {j: < j < p- 1, mj / 0}. Then 21' := H}^U {fj} jE n 2 is 
an orthonormal basis in H. Set 2tg := {fj}jen 2 - Observe that card(2l(j) = 5. 
The k-th element, counting from zero, of the set 2lg , arranged in the order 
of construction of its elements, we denote by v T+k , k = 0,1, ...,5 — 1. Then 

a , = {«*}a" 1 .^ = K}S" 1 - 

Denote by A4i^(&) the matrix of the operator Eh — C{A © in the 
basis 21, £ £ B. Here is an analytic in ID operator-valued function which 
values are linear contractions from H © D{A) into H R(A). Then 

^i,cw = (a^-c(^©^c)]^%y;S 1 = ( ^ ^ ), 

where 

A),c = {([ E H-C(A®^ c )]u k ,u j ) H ) T j ~ =0 = {[uk-C,Au k ,Uj) H y~ k={) 

= Ir-C{{Vk,U j ) H y~ k l =Q , (15) 

B 0X ($) = (([E H -C(A®$ c )}u k ,u j ) H ) ^ T _ lj 

= {{U k - C^C^'^)if) <i<r-l, r<fc<r+5-l 
= _ C ((^C U fc>%)i27 0<i<T-l, r<fc<r+5-l ' 

Co, c = {{[E H ~ C(A © $ c )] Ufc , Ui ) ff ) T < 

= (K - c^«fc,^) if ) r < j < T+5 _ lj 

= -C{(vk,Uj) H ) T ^ T+s _ h o^^!, (16) 

L>0 )C ($) = {([EH-aA^^^^^r^r+S-l, r<k<r+6-l 
= ((U k - C^C^'^)/f) T < i < T +5-l, r<k<r+S-l 

= I S -C (($C u fc> u j)ff) T <j< r +,s-i, T < fc < T+i _! , C € B. 
Observe that the matrix is invertible, since it is the matrix of the 
operator P D ( A ){Eh - (A)P D ( A) = E D( ^ A) - (P D (A)AP D ( A ), considered in 
the Hilbert space D(A), with respect to 2I2, C 6 B. Notice that matrices 
A),c>Co,c> C 6 B, can be calculated explicitly using relations © and flU). 

Denote by £ £ B, the matrix of the operator $£, acting from 
-0(A) into -ff © R(A), with respect to the bases 2I3 and 21' 3 : 



Ft 



ii • 
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Then 

T+S-l 

$(Uk= Yl h( l ' k ) v ii r<k<T + 5-l, 

l=T 

and 



/ /t+6-1 n 



H' 0<j<T-l, r<k<T+5-l 



/t+S-1 \ 

= "C £ ( W ) H / c (J,fc) , CGB- 

V J=T / 0<j<T-l, T<fc<T + <5-l 

Set 

W := ((^,^) H )o<i<r-l, r<Kr+5-l " ( 17 ) 

Then 

B 0t( ($) = -CWF C , C G 

We may write 

d , c ($) = j (5 -c(( £ /c(i^K^) j 

\ \ l=T ' H/ T<j<T+S-l, T<k<T + 8-l 

/t+S-1 \ 

= is-C[ £ (^iWc(^)J > Cei- 

V i=T / T<j<T + S-l, T<k<T+8-l 

Set 

T := ((11,11;)^^ T <t< T+S -i • ( 18 ) 

Then 

D 0>f (*) = /* - CTF C , (6»- 

Thus, we may write 

where A),o Co.C are given by (USD,®, and TV, T are given by ([T7|) .([TS ]) . 
Let apply the Probenius formula for the inverse of a block matrix [H p. 59]. 
Then 

m-u*) = ( ~ c^flTH*)^ * 



s 



±A+ C -<A+ C WF C H^(<S>)C , C A+, * 



C 
* 



where by stars (*) we denoted the blocks which are not of interest for us, 
and 

H C ($) = h- (TF C + CCo, C A^WF c = I 5 - (TF C + ^C , c i+ ( W ( 



= h + (£-C , c A+ c W - (T^j F C , (eO. 



(19) 



Here Aq^ denotes the adjoint matrix of i.e. the transpose of the co- 
factor matrix, and 

h c = detA), c . (20) 
Let (£D. The minor of A4^^(&), standing in the first p rows and the first 
p columns, we denote by M.2,c(^)- The minor of Aq^, standing in the first 
p rows and the first p columns, we denote by A\^. The first p rows of Aq^ 
we denote by A 2 ^- The first p columns of Aq^ we denote by ^3,^. Then 

M 2 , C ($) = ^A U -^A U WF C H~H$)C , C A U , (21) 

Observe that M2,((®) is the matrix of the operator P Ln [E h - ((A © <£>f)] _1 Pl n , 
considered as an operator in Ljv, with respect to the basis 2li, ( £ D. 
Consider the following operator from to L^: 

N-l N-l 

K y c n e n = ^ ^ c n x n , c n £ C, 

n=0 n=0 

where e n = (6 nt o,Sn,ii ■ ■ ■ ,8n,N-i) £ C N . Let K, be the matrix of K with 
respect to the orthonormal bases {e n }^~Q and Sti: 

/C = ((Kek,Uj) H ) Q <j< p _^ < fc < A r_ 1 = ((Xk, U 'j) h) 0<j<p-l, 0<k<N-l ■ ( 22 ) 

Then we may write 

= (p Ln [E H - aA®^)]- 1 P LN Ke k ,Ke 3 ) h 
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K*P Ln [E h - ((A ^c)]" 1 P LN Ke k , e/)^ , ( G 1. 

Observe that the right-hand side is equal to the element of the matrix 
/C* M. 2 ,£(3>)/C , standing in row j, column k. By ([8]) we may write: 



2n -1 

—— - t dM T {t)=K,*M u (<$>)1C, CGB. (23) 
Set 

C c = £C 0X A+ C W - Ch c T, A c = /CM U /C, 
B C = JC*A U W, B c = C 0iC A 3X JC, (eB. (24) 
By (USD,®,® we get 

r t^" mt v = ^ - ^ ( /s + ir c c ^) ' d <- 

where E D. 

Theorem 2 Let t/ie moment problem (CP, u>i£/i d > 1, be given and condi- 
tion with Td from (d|), be satisfied. Suppose that the moment problem 
is indeterminate. All solutions of the moment problem |7]) can be obtained 
from the following relation: 



2tt j 



1-Ce 



dM T (t) 



= /^ Ac "^ B ^( /5+ /^ C?i?C ) ° C ' CGD ' (25) 

where A^, B^, C^, D^, are matrix polynomials defined by \21$ , with val- 
ues in CnxN, Qvx<5; CsxS, CsxN, respectively. The scalar polynomial h^, 
deg/i^ < t, is given by Ii20\) . Here Fs is an analytic in B, CsxS-valued func- 
tion which values are such that F^F^ < 1, V£ € B. Conversely, each analytic 
in B ; C$xs-valued function such that F£F^ < 1, VC £ B ; generates by rela- 
tion [25\) a solution of the moment problem |7J). The correspondence between 
all analytic in B, C$ x s-valued functions such that F£F^ < 1, V£ £ B, and 
all solutions of the moment problem (CP is bijective. 

Proof. The proof follows from the preceeding considerations. □ 
Example 2.1. Let N = 3, d = 1, S 
Consider the TMTMP with moments Sq,S\. It is straightforward to check 
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that condition ([3]) holds and condition (C) of Theorem Q] fails. Thus, the 
TMTMP is solvable and indeterminate. The matrix T\ from ([2]) has the 
following form: 

1 10 1 1 \ 
110 110 
1 
110 110 ' 
110 110 
1 / 

Let H be the Hilbert space described above, after formula ([5]), and {x n } n _ 
be elements with the property ([5]). 

Let us apply the orthogonalization procedure (fTTj) . (fl~2l) to the elements 

X ,Xl,X2,X 3 ,X A ,X 5 . 

Step 0. Calculate no = ||#o||i? = y ( x o, %q)h = \/7o,o = 1/0. Then we 
set y = Xq. 
Step 1. Calculate 

n \ = ||«i-(»i>2/o)fl-yo||H = lki-0ci,^o)i^o||H = (»i -71,0^0) a?i-7i,oa?o)fr 
= (x 1 - x ,xi - x ) H = (x 1 ,x 1 ) H - (xi,x ) H - Oo, x\)h + (%o, x )h 
= 7i,i - 7i,o - 7o,i + 7o,o = 0. 

Therefore we pass to the next step. 
Step 2. We calculate 

n 2 = 11^2 - (x2,yo)Hyo\\H = \\X2 ~ 72,0^o||/f = (x 2 ,X 2 ) H = 72,2 = 1- 

Set 

y 2 = x 2 - (x 2 ,yo)Hyo = x 2 - 72,02/0 = x 2 . 

In step 3 we obtain 77,3 = 0, in step 4 we get 714 = 0. Finally, in step 5 we 
get n 5 = 1, and 2/5 = x 5 . 

Set 21 = {2/0,2/2,2/5}- Let u := y = x , ui := y 2 = x 2 , u 2 := y 5 = x 5 . Then 
21 = {uk}\ = Q- Observe that in our case we have: p = r = 2, S = 1. 
Set t>o := ^4uo = ^4xo = £3, v\ := Au\ = Ax 2 = x$. Let us apply the 
orthogonalization procedure (fl~3j) . (fill) to the elements Vq, v%, Uq, U\. 
Step 0. Calculate 

m l = \\uq-{uq,vq) h vq-{uq, v^hViWh = \\xo-(x ,x 3 ) h x 3 -(x ,x 5 )hx 5 \\ 2 h = 

\\XQ - 10,3X3 - 10,5X 5 \\h = \\XQ - X 3 \\h = (XO - X 3 ,X - X 3 ) H 



(7, 



n,m) n ,m=0 
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= 7o,o - 70,3 - 73,0 + 73,3 = 0. 

Then we pass to the next step. 
Step 1. Calculate 

m\ = \\ui - (u 1 ,v )hv - (iti,ui)fl-«i[|H = 

\\x 2 - (x 2 ,X 3 ) H X 3 ~ (x 2 , X 5 ) H X 5 \\h = \\X 2 ~ 72,3^3 ~ 72,5«5||# 
= II^Hh = ( x 2, X 2 )h = 72,2 = 1. 

Set 

h = u\- (ui,v )hv - (ui,vi) H vi 

= X 2 - (x 2 , X 3 ) H X 3 - (x 2 ,X 5 ) H X 5 = X 2 - 72,3^3 ~ 72,5^5 = X 2 . 

Set 21' = {v ,vi, fi}. Let v 2 = fi = x 2 . Then 21' = {v k } 2 k=0 - 
By (H7J),(1IED we may write 



(x 2 ,x ) H \ _ ( 72,0 \ _ ( 
72,2 / V 1 



T = ({vi,Uj) H ) 2 ^.^ = (v 2 ,u 2 ) H = (x 2 ,x 5 ) H = 72,5 = 0. 

By (US]),® we calculate: 

i-C o 
o 1 



4>,c = h 



C ,C = ~C {(vki u j) H ) 2 <j< 2 , o<k<i = -C(73,5,75,s) = -C(0,1). 

Then h c = det A 0>c = 1 - C, >4+ f = ( J 1 - C ) = = ^ 2,C = 
By (I22p we may write 



7o,o 7i,o 72,o \ _ / 1 1 



£ (( z fc>%');j)o<i<i, o<fc<2 «, - n2 -, | 2 v, 2 y y (, o j 
Using (|24l) we calculate 

1 1 \ / 

C C = -C 2 (1-C), A c = ( 1 1 , B c = 

o o i-c J V i-C 
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-C(o,o,i-c). 



Finally, by (|25|) we obtain: 



2t 



: dM T (t) 



o 1 - Ce f 



/ i-c i-C 

i-C i-C u 

V o o i + C 2 t4% 



In particular, if we choose Fs = 1, then 

/ m(i) m(i) 

M(i) = m(t) m(i) 

\ m(t) 

where 



m(t) 



0. 



if t = 



1, ifte(0,27r] 
is a solution of the TMTMP. 



, m(t) 



t G [0, 2vr] 



C G D. 



0, 



if t = 



|, ift€(0,7r] 
1, iftG(7r,27r] 
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This paper is a continuation of our previous investigation on the trun- 
cated matrix trigonometric moment problem in Ukrainian Math. J., 2011, 
63, no. 6, 786-797. In the present paper we obtain a Nevanlinna-type formula 
for this moment problem in a general case. We only assume that we have 
more than one moment, the moment problem is solvable and the problem 
has more than one solution. The coefficients of the corresponding matrix 
linear fractional transformation are explicitly expressed by the prescribed 
moments. Easy conditions for the determinacy of the moment problem are 
given. 
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